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Phase transitions in systems described by Bose-Fermi-Hubbard model on a lattice with two nonequivalent 
sublattices are investigated in this work. The case of hard-core bosons is considered and pseudospin for- 
malism is used. Phase diagrams are built in the plain of chemical potential of the bosons-bosonic hopping 
parameter. It is shown that in the case of anisotropic hopping, the region of the supersolid phase existence is 
possible for a smaller parameter space. 
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1. Introduction 

Theoretically, the systems of ultracold atoms are well described by the Hubbard model. For 
the case of bosonic atoms, the investigations of the Bose-Hubbard model predicted a superfluid to 
Mott insulator transition in ultracold atomic gases |T1[1|- The experimental verification of this phase 
transition was done by Greiner et al. [3J and the study of this model has rapidly developed since 
then. Besides the bosonic atoms, the mixture of bose-fermi atoms in optical lattices was realized 
and the Bose-Fermi-Hubbard model (BFHM) is under intensive investigations f^-'lO]. As concerns 
the BFHM, phase diagrams are more complicated because the effective interaction between bosons 
is generated due to the presence of fermions and this leads, for example, to the possible appearance 
of a supersolid (SS) phase (which is characterized by a simultaneous presence of a density wave 
and phase order in the condensate). 

The lattice depth, dimensionality, geometry, and filling factor can all be controlled. Apart 
from the uniform lattice potentials, other lattice topologies can be realized. For example, the 
superposition of two standing-wave lattices of different wavelengths leads to a superlattice with a 
spatial modulation of the lattice well depths. In pTj an optical lattice of double wells was realized 
by combining laser beams with the in-plane and out-of-plane polarizations of light. The properties 
of the system of bosonic atoms in superlattice potential were studied in |12ffT7] , where the presence 
of the modulated potential leads to the spatial modulation of the bosonic concentration. In [TB] the 
zero temperature phase diagram of binary boson- fermion mixtures in one-dimensional superlattices 
was investigated using an exact numerical diagonalization technique. It should be noted that almost 
all the mentioned studies were restricted to the case of zero temperature. In our previous work fTU] 
we analyzed phase transitions in the BFHM at finite temperature. 

The Bose-Fermi-Hubbard-type model can also be applied in the description of intercalation of 
ions in crystals (for example, lithium intercalation in Ti02 crystals). Ion-electron interaction can 
play a significant role in these systems. At intercalation, the chemical potential in such systems 
is displaced into the conduction band. Phase separation into Li-poor and Li-rich phases occurs 
in such crystals and this two-phase behaviour leads to a constant value of the electrochemical 
potential [191 [20] (this is taken in consideration when constructing batteries). 
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2. Model and results 



In this work we consider the thermodynamical properties of the BFHM on a lattice with two 
nonequivalent sublattices. We consider the hard-core limit (infinite on-site boson-boson interac- 
tion) . The geometry of the considered lattice is shown in figure [l] where we take into account two 
parameters of the fermionic hopping (ti and and bosonic hopping {Tti and fl2)- This type of 
lattices can be experimentally realized. For example, the lattice with similar topology was experi- 
mentally realized in 



t2 




(a) (b) 

Figure 1. (a) Schematic picture of a double-well lattice in 3D case, thin lines with arrows 
correspond to hopping with parameter t2 and thick lines correspond to hopping with parameter 
ii. (b) Double- well potential corresponding to the cross section for the dashed line in (a), x 
direction. 



Using the pseudospin formalism, the Hamiltonian of the model is written in the following form 



(2.1) 



The pseudospin variable Sf = 1/2 when a boson is present in a site i (boson concentration rtb = 
5*^ + 1/2), cf and are fermionic creation and annihilation operators, respectively. The first 



and the second terms in equation (2.1 1 ) are responsible for nearest neighbour boson and fermion 
hopping, respectively; g-term accounts for the boson-fermion interaction energy. We consider the 
grand canonical system and introduce the bosonic and fermionic chemical potentials h and fj,, 
respectively, a — A, B is sublattice index. Notice that — h, fiA = fJ-, h^ — h — E.fiB^fJ'^E 
(see figure [T]) . 

We use the mean field approximation (MFA): 

ns+sj ^ n{s+)sj + ns+{sj)-n{s+){s-), (2.2) 

which is appropriate in the case of weak boson-fermion interaction when the fermionic band is not 
split due to this interaction. Application of the mean field approximation to strongly correlated 
systems in the limit of a weak on-site correlation makes it possible to satisfactorily describe their 
properties. For example, application of the similar approximation to the fermionic Hubbard model 
in the limit of the weak on-site electron correlation makes it possible to describe its magnetic 
properties. In our previous papers (see, for example, jH]), we considered the case J7 = (in this 
case the model is similar to the Falicov-Kimball model) and decoupled the on-site interaction term 
in a similar way. We showed that in the case of weak coupling, our results qualitatively agree 
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with those obtained within the framework of the dynamical mean field theory for the Falicov- 
Kimball model. In addition, in the case of the Bose-Hubbard model, the kinetic energy term is 
often considered within the mean field approach. 
The Hamiltonian, therefore, is written as follows: 



H = -2nY^{St^{SI,)+Sl{Sl))+N^l{Sl){Sl)~Y.^,,,^^ 



Here, N is the number of lattice sites and O — rjj^jg (we can see that in the applied approxi- 
mation the difference between fli and Q2 does not play any role). 

At first we diagonalize the pseudospin part of the Hamiltonian using the unitary transformation 
in the pseudospin subspace: 



Sa = COS 9a + CTq sin 9a , 

S'a = (TlcOs9a~ a'aSlli 9a., 

. . 2fi(5|) ha- gua 

Sm 9a = r , cos 9a = r 

An. A„ 



Xa = ^{gna-haY + {2n{Sfi)\ 



ET5(5^)K>+iVf^(5X)(^i), a,p^A,B, (2.3) 



here we passed to k- representation, tk = |ik|c''^. In one-dimensional case |ikP = ^1+^2 + 
2tit2 cos 2k, cos 7 = ■^^^-'-j^^p^ , sin 7 = ^^p^- Similarly, in three-dimensional case for the con- 
sidered lattice geometry (see figurejlj |tkp = ti+t2 + 2tit2Cos2kx+4:t2{cosky + coskz)coskx{ti + 
^2) + 4i2(cosfcj, + cosfc^)^. To diagonalize the fermionic part of the Hamiltonian we perform the 
following unitary transformation 

CkA = flkA cos -I- OkB sin (j)e^'^ , 

CkB = — akA sin 0e~'^ -I- OkB cos (2.4) 
cos 2^ . /^g-^^A + i^ ,.^20 ^ -l^.l (2.5) 



It should be noted that fermionic spectrum is always split when ti 7^ t2, this is true even in 
the case when {S^) — (S^) (when (S^) 7^ (S^) there is an additional splitting with the gap 
sK'^a) ~ ('S'b)I llQi)- one-dimensional case, this splitting is of the order of 2\ti — t2\- In what 
follows we will discuss the effect of the splitting on the phase transition picture. 

Now we can write 

H = E^k.aJ,„ak,a-E^"^--Efff(^")("")+^"<^A)(^§>, 

. + (58) + E ^ - .(Si) - g )' ~, 
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Now we are ready to obtain the system of equations for bosonic and ferniionic concentrations 



, , 1 I l+cos(2(/)) 
= mI^S 2 



ha - , ( (3Xa 

_^tanhl — 




exp 



Xk/3 — fJ' 

T 



-I -1 



(2,7) 



The grand canonical potential can be written as [TIJ] 



7V72 



T 



El- 



1 + exp 



A 



kA 



-Tin ( 4 cosh ^^—^ cosh 



T 



1 + exp 



-5((nA)(5i) + (nB)(^^)) + 217(51) (2.8) 



It was shown in [10 that coming from the set of the mean-field equations (2.7) we can obtain 
the condition of the appearance of nonzero values of Sn = {ua) — ("-b); SS'^ = (5*1) — (S^), 
SS'^ = (5*1) — (Sq) (which play the role of the order parameter for the modulated phase) and this 
condition coincides with the condition when the static density-density correlator {S^ S^)q=Tri^^ = 0) 
calculated in the random phase approximation diverges. 



In what follows we consider a three-dimensional case. Coming from the set of equations (2.7) 



it can be shown that (S^) ~ when $7 < 2T. Therefore, at finite temperature we can consider the 
transition from the uniform nonsuperfluid phase (at low temperatures this is Mott insulating (MI) 
phase) to the charge density wave (CDW) phase for small values of the bosonic hopping parameter 
(i7 < 2T). We use the equations for averages (2.7) and the expression for the grand canonical 
potential (|! 

J7) and the dependencies of 



to find thermodynamically stable states. 
2] |3] [4]we show the phase diagrams in the plane (h 



In figures 

the bosonic and ferniionic concentrations on bosonic chemical potential at low temperature for 
different cases when ti = t2, ti ^ t2, E = 0, E 0. We consider the regime of the fixed fermionic 
chemical potential and from figures [2] [3] [4] we can see that the phase transition from the uniform 
to chess-board phase can be of the second (solid line) or first (dashed line) order. The existence 
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Figure 2. (a) Phase diagrams in the {h — D.) plane for ti = t2 = 1/4:, g = —0.4, E = 0, 
jj, — 0, T — 0.005. Solid (dashed) lines denote the second (first) order phase transition lines, 
(b) Dependencies of the bosonic and fermionic concentrations on the chemical potential of bosons 
at n = 0.2. 
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Figure 3. (a) Phase diagram in the {h — f2) plane for t\ = 1/4, t2 = 1/6, g = —0.4, E = 0, 
jj, — 0, T = 0.005. Solid (dashed) lines denote the second (first) order phase transition lines, 
(b) Dependencies of the bosonic and fermionic concentrations on the chemical potential of bosons 
at n = 0.21. 



of the phase transition of the first order leads to phase separation (in the regime of the fixed 
concentrations) into the uniform and CDW phases (this was illustrated in [5^). The presence of 
anisotropic hopping (ti ^ t-i) leads to the narrowing of the region of the supersolid (SS) phase 
(this phase is characterized by the simultaneous presence of a density wave and phase order in the 
condensate) . This is due to the appearance of the above mentioned gap in the fermionic spectrum 
leading to the plato-like behaviour in the dependence of the fermionic concentration on the bosonic 
chemical potential in the SF phase, see figure[3](b). In our numerical calculations we considered the 
case <i = 1/4, ^2 = 1/6, but similar conclusion (about the narrowing of the SS region) is also valid 
for other values of ti, t2 (for example, we also calculated the phase diagrams for the cases t\ = 1/3, 




(a) (b) 

Figure 4. (a) Phase diagram in the {h — Vt) plane for tx = 1/4, t2 = 1/4, g = —0.4, E = 0.1, 
/i = 0, T = 0.005. Solid (dashed) lines denote the second (first) order phase transition lines, 
(b) Dependencies of the bosonic and fermionic concentrations on the chemical potential of bosons 
at n = 0.2. 
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^2 = 1/4, and ti = 1/6, t2 — 1/4, but here we do not present these diagrams because they are 
similar to that shown in figure [3|. With an increase of temperature, the first order phase transition 
transforms into the second one and then disappears at some critical value of the temperature, as it 
was discussed in our previous work [TU] at the construction of T — h diagrams in the case ti ~t2- 

The case E illustrated in figure |4j As 

we can see, the presence of the modulated po- 
tential leads to the modulation of bosonic and 
fermionic concentrations in sublattices (at any 
finite value of E we always have such a modu- 
lation), we call this phase a normal (NR) non- 
superfiuid phase to distinguish it from a true 
charge-density- wave phase, for which the trans- 
lational invariance is spontaneously broken due 
to the presence of boson-fermion interaction. 
In our approach, we can distinguish between 
these two phases (NR and CDW) only when 
the phase transition between them is of the first 
order with jumps of the fermionic and bosonic 
concentrations. As we can see, the phase dia- 
gram is asymmetric in the case E' 7^ in com- 
parison with the case E = Q (see figures [2] |3| 
because the chemical potential of particles on 
one of the sublattices is shifted away from the 
half-filling case. When we increase the value of 
the double-well potential E, the line of the first 

order phase transition between modulated phases at higher values of the bosonic chemical potential 
{h « —0.03 at il = in figure |4] (a)) gets shorter and the line of the first order phase transition 
at lower values of the bosonic chemical potential {h w —0.25 at = in figure [4] (a)) gets longer 
(here we do not present phase diagrams for other values of E because this tendency is also seen in 
figure |4] (a)). 

The presence of the first order phase transition is connected with the effective interaction 
between bosons via fermions. As the temperature increases, this effective interaction decreases and 
at some critical value of temperature this first order phase transition disappears (see figure [5]) . 
At low temperature the bosonic concentration is almost the same in both sublattices (nb = or 
TT-b = 1 depending on the value of the bosonic chemical potential) when we go into NR phase (at 
zero temperature in the case of the Bose-Hubbard model this phase is sometimes referred to as a 
hole vacuum or particle vacuum [14 when nb = 1 or nb = 0, respectively), see figure |4](b). We 
do not use the notation "SS phase" in figure |4] because we cannot distinguish between SS and SF 
phase (for the case of a true SS phase, the translational invariance is spontaneously broken rather 
than due to the presence of the modulated potential). 



Figurs 5. Phase diagram in the (T — h) plane 
for tx = 1/4, t2 = l/4:,g = -0.4, E = 0.1, a* = 0, 
Q, = 0.2. Solid (dashed) lines denote the second 
(first) order phase transition lines. 



3. Conclusions 

The Bose-Fermi-Hubbard model on a lattice with two nonequivalent sublattices has been in- 
vestigated in this work, and this type of optical lattices can be realized experimentally (see, for 
example, [H]). We considered the hard-core boson limit and investigated the thermodynamics of 
the model for the case of weak boson-fermion interaction. The phase diagrams in the plane of the 
chemical potential of bosons-bosonic hopping parameter were built and it was shown that the tran- 
sition from the uniform to modulated phase can be of the first or the second order (in the regime 
of the fixed concentrations this leads to phase separation). The gap in the fermionic spectrum 
appears due to the presence of anisotropic hopping and this causes the narrowing of the parameter 
region of supersolid phase existence. The presence of the effective interaction between bosons via 
fermions leads to the first order phase transition between modulated phases at the presence of the 
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double-well potential E. It should be noted that we considered thermodynamics in the mean field 
approximation and our future investigations will be devoted to the corrections beyond the mean 
field approximation. 
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Mofle/ib Bo3e-0epMi-ra66apAa hb r'paTi^i is abomb 
HeeKBieajieHTHMMM ni/^r'paTKaMM 

T.C. MMCaKOBMH 

iHCTi/iTyT 4>i3i/iKi/i KOH/^eHCOBaHi/ix CMCreM HAH VKpaiHi/i, Byji. I. CBeHL^iLfbKoro, 1, 79011 JlbsiB, VKpaiHa 

y i^iii po6oTi /^ocniflweHO 4>a30Bi nepexoAi/i y CHCreMax, u^o onwcyiOTbOi MOfl,eAj\K> Bo3e-OepMi-ra66apfla 
Ha rpaTLfi 3 pfioma HeeKBisajieHTHMMH ni^rpaTKaMi/i. Po3rji;iHyTO Bi/inaflOK >KopcTKi/ix 6o30hIb ra Bi/iKopi/i- 
craHO nceBflOcniHOBi/iu c|3opMajii3M. FloGyAOBaHO c|}a30Bi /^iarpaMH y njioii^HHi xiMiHHi/iM noreHL^iaji 6o30HiB- 
6030HH1/1M napaMerp nepecKOKy. noKa3aHO, mo y BiinaflKy aHi30TponHoro nepecKOKy o6yiacTb icHysaHH;) raK 
3BaHoY "cynepcojiifl" 4>a3H 3By)KyeTbC5i. 

K/hohobI c/iOBa: MOflejib Bo3e-0epMi-ra66apfta, o/tthwh; tparKi/i, (pasoBi nepexoflu 



43301-8 



